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We report on the new development in studying 
supersymmetric quantum mechanical systems. 
The main goal is to solve the ten (space-time) 
dimensional (D=10) supersymmetric Yang-Mills 
quantum mechanics (SYMQM) with the SU(N) 
gauge group [Q. At infinite N it provides rel- 
atively simple model of M-theory|^][]. In spite 
of recent results that question the exact equiva- 
lence between the two, D=10 SYMQM possesses 
lot of fascinating properties required in M-theory 
(e.g. threshold bound states). Independently, 
these systems provide a good laboratory to study 
various consequences of supersymmetry. 

We use the hamiltonian formalism adapted to 
supersymmetric systems with local gauge invari- 
ance j|. It consists of two steps: 1) creation of a 
finite basis of gauge invariant states, and 2) alge- 
braic calculation of the matrix representation of 
the hamiltonian, followed by its numerical diago- 
nalization. Algebraic calculations are automated 
by implementing standard rules of quantum me- 
chanics in an symbolic language like Mathemat- 
ica. Faster, compiler based, version is now also 
available. Given the spectrum and energy eigen- 
states, any other relevant observable can be read- 
ily obtained using our computer based rules of 
quantum mechanics. 

Such an approach necessarily introduces a cut- 
off - the size of the basis. We choose for it the 
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3 For finite N it represents compactified version of M- 
theory. 



gauge and rotationally invariant total number of 
bosonic quanta B — a l a a^ . Hence the cut off 
basis consists of all states with B < N cut and 
all allowed fermions. Changing the size of the 
basis gives the quantitative measure of the finite 
./Vcut effects similarly to lattice calculations. 

The technique has been applied to Wess- 
Zumino quantum mechanics, D=2 and D=4 
SYMQM and to D = 5 - 10 YMQM, all based 
on the SU(2) gauge group ||. In all cases studied 
until now the spectrum of lower states converges 
with -/Vcut well within the reach of present com- 
puters. Clear SUSY signatures are seen. Simi- 
lar methods have been independently developed 
in Q for lower dimensional supersymmetric field 
theories^ Since || a notable progress, in increas- 
ing the Fock space, has been achieved by M. Cam- 
postrini. Here we report the new results for four 
(space-time) dimensional quantum mechanics. 
D=4 SYMQM. 

The system is described by nine bosonic coordi- 
nates x l a {t), i = 1, 2, 3; a = 1, 2, 3 and six indepen- 
dent fermionic ones contained in the Majorana 
spinor ip"(t), a = 1, ...,4. Hamiltonian reads 

1 - - q 2 .... 

= 2 PaPa + -^ e *bc£adeX b X J c X d X J e (1) 

with the Dirac a matrices. We work in the Ma- 
jorana representation of Itzykson and Zuber. 

The system has rotational symmetry, 
with the spin(3) angular momentum J 1 = 

4 See Uwe Trittmann contribution. 
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the gauge invariance with the generators 
G a = fabc (x^Pc ~ f^bVc) > and supersymme- 
try generated by the charges Q a — {T k ^p a ) a p k l + 

Generation of basis, construction of Majorana 
fermions and details of the calculation of hamilto- 
nian matrix and other observables are described 
in H . Here we present new results for higher cut- 
offs ( N cut < 8 ) in all fermionic sectors. Hamil- 
tonian (||) conserves fermion number F = f^fai 
hence it suffices to diagonalize H separately in 
each fermionic sector. Moreover, as a conse- 
quence of a particle-hole symmetry there are only 
four independent sectors. 
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Figure 1. Cut-off dependence of the lower levels 
in all four fermionic sectors. 



F = 1 we again observe fast convergence with 
increasing 7V cu t - our bases are adequate to ex- 
tract the "infinite volume" limit of lower glue- 
balls and gluinoballs in both sectors. One also 
observes emergence of the supersymmetric multi- 
plets spanned across different fermionic sectors^. 
In "fermion rich" sectors F — 2,3 situation is 
even more interesting. In addition to quickly con- 
vergent with N cut states we observe levels with 
strong iV cu t dependence whose energies tend to 
zero. In fact the cutoff dependence can be used 
to identify discrete and continuous spectra. It was 
found that localized, discrete bound states have 
fast (e.g. exponential) convergence with iV cu t , 
while continuous, non localized spectrum is sig- 
naled by slow 0(1/N cut ) dependence. Therefore 
our results confirm (and quantify) famous SUSY 
predictions of the continuous spectrum in this 
model We also establish that the SUSY vac- 
uum is the lowest, J = 0, state in the F = 2 
sector, that it belongs to the continuum hence is 
non normalizable, and in fact is doubly degener- 
ate (F = 2 and 4 sectors) . 
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The spectrum of the theory (Fig.l) is rather 
rich and behaves differently in different sectors. 
In F = sector our spectrum converges to 
the well known classic results of Liischer and 
Munster for zero volume glueballs || providing 
satisfactory test of the whole approach^]. For 

5 Together with van Baal we have also found complete 
agreement (for F = and F = 2) with the recent adapta- 
tion of his method of supersymmetric spherical harmonics 

I- 



Figure 2. Witten index of D=4 SYMQM for up 
to 8 bosonic quanta. 



Witten index I W (T) = S^-l)^ exp {-TEi), 
provides a global measure of the restora- 

6 SUSY multiplets can be identified in detail using super- 
symmetric charges and angular momentum operator. 
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tion/violation of SUSY. It is discontinuous at 
T = with lw(0) depending on the regulariza- 
tion. For discrete spectrum it is time independent 
and equals to the number of ground states. For 
continuous spectrum it may have a mild T depen- 
dence, still Iw(oo) gives the number of localized 
ground (unbalanced by SUSY) states. The in- 
dex shown in Fig. 2 was obtained directly from 
our spectrum. Compared to || we observe fur- 
ther slow approach to the large iV cu t limit, con- 
sistent with all above properties. The physical 
(bulk) value lim^^o Iw (T) seems to point to- 
wards 1/4 calculated from nonabelian integrals 
H and marked by the flat line. There is no indi- 
cation of the threshold bound state in agreement 
with Fig.l and with theoretical expectations. 
From jour to ten space-time dimensions. To as- 
sess the feasibility of our approach to the full 
D=10 theory we have calculated the spectra of 
the purely bosonic hamiltonians in all interme- 
diate dimensions 4 < D < 10. Thereby results 
shown in Fig. 3 extend for the first time the spec- 
trum of the zero volume glueballs to higher di- 
mensions. This generalization is straightforward 
in our approach, however it demands larger, but 
realistic, computer effort. In all dimensions we 
see fast convergence with iV cu t indicating the lo- 
calized nature of glueball bound states. Since 
our cutoff respects rotational symmetry for all 
d=D-l, we observe exact SO(d) degeneracies in 
the spectra. The ordering of the first three mul- 
tiplets, shown in Fig. 3, remains the same as in 
d = 3: scalar, two dimensional symmetric tensor 
with dimension g — d(d + l)/2 — 1, and scalar 
again. The largest uncertainty (i.e. that of the 
third d=9 level) is around 25% at present and di- 
minishes gradually to below 1% at d = 3. Turn- 
ing on d as another parameter allows for further 
confrontation with analytical calculations. For 
example the simple eye-ball estimate from Fig. 3 
suggest the power dependence E ~ d~ 15 to be 
compared with the mean field result at low tem- 
perature d 4 / 3 J|]. 

We conclude that calculations in higher dimen- 
sions including d — 9 are perfectly realistic even 
if more time consuming. Consequently a path 
towards the quantitative study of the full D=10 
supersymmetric theory is now available. 
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Figure 3. Zero volume glueballs in higher dimen- 
sions. 
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